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PREFACE 


This  paper  describes  the  detection  coil  design  and 
construction  involved  in  making  an  earth's  field  magneto- 
meter which  has  substantially  a  continuous  information 
rate.   Background  information  concerning  the  proton 
free  precession  magnetometer  is  supplied  and  the  test 
results  of  the  high  information  rate  system  are  presented 

The  authur  worked  on  this  project  in  conjunction 
with  Lieutenant  Grover  M.  Yowell  at  Varian  Associates, 
Palo  Alto,  California,  in  the  early  months  of  1956. 

The  author  wishes  to  thank  Varian  Associates, 
Dr.  Martin  Packard,  and  all  the  members  of  the  Instrument 
Division  whose  cooperation  made  the  accomplishment  of 
this  project  possible. 


ii 


TABLE  OF  CONTENTS 


Item 


Title 


Page 


Chapter  I 
Chapter  II 
Chapter  III 
Chapter  IV 
Chapter  V 
Chapter  VI 
Chapter  VII 
Chapter  VIII 
Bibliography 


Introduction 1 

Nuclear  Magnetic  Resonance  ....  6 

Free  Precession  Magnetometer  ...  17 

Coil  Design  and  Construction  ...  28 

Coil  Field  Measurements 45 

System  Design 58 

System  Test 61 

Conclusion 69 

72 


iii 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

2-1     Proton  Alignment  in  Steady  Field  9 

2-2     Proton  Precession  in  Suddenly  Applied  Field  9 

2-3     Nuclear  Magnetic  Resonance  Apparatus  10 

2-4     "v"  mode  Signal  12 

2-5     "u"  mode  Signal  12 

2-6     Portion  of  Ethyl  Aloohol  Spectrum  12 

3-1     Free  Precession  Signal  20 

3-2     Station  Magnetometer  Block  Diagram  21 

3-3     Station  Magnetometer  22 

3-4     Polarizing  Coils  and  Sample  24 

3-5     Station  Magnetometer  Record  27 

4-1     Magnetic  Field  from  a  Straight  Wire  32 

4-2     Magnetic  Field  from  Two  Straight  Wires  32 

4-3     Components  of  Double  Solenoidal  Coil  38 

4-4     Assembled  Double  Solenoid  39 

4-5     Toroidal  Coil  43 

5-1     Field  Measurements  Using  Station  Magnetometer  46 

5-2     10,000  Turn  Pick  Up  Coil  47 

5-3     External  Magnetic  Field  Intensity  of  Double 

Solenoid  1       ,  50 

5-4     External  Magnetic  Field  Intensity  of  Double 

Solenoid  2  51 

5-5     External  Magnetic  Field  Intensity  of  Double 

Solenoid  3  52 

iv 


>» 


LIST  OF  ILLUSTRATIONS 
Continued 


Figure  Page 

5-6     External  Magnetic  Field  Intensity  of  Toroidal 

Coil  1  53 

5-7     External  Magnetic  Field  Intensity  of  Toroidal 

Coil  2  54 

5-8     External  Magnetic  Field  Intensity  of  Toroidal 

Coil  3  55 

6-1     Dual  Coil  System  Simplified  Block  Diagram  59 

7-1     Double  Solenoid  Arrangement  for  System  Test  62 

7-2     Noise  Induced  Fluctuations  63 

7-3     Data  Presentation  1  66 

7-4     Data  Presentation  2  67 

7-5     Data  Presentation  3  .  68 


CHAPTER  I 
INTRODUCTION 

The  measurement  of  the  earth's  magnetic  field  has  been 
of  interest  to  scientists  and  navigators  for  many  years. 
More  recently  such  data  has  become  useful  to  geophysicists 
and  to  the  military.  Field  measurements  have  been  made 
with  compass  needles,  earth  inductors  (1),  saturable  core 
devices  (2),  (3),  electron  beam  deflection  tubes  (4), 
electron  paramagnetic  resonance  (5),  and  by  free  precession 
nuclear  magnetometers  (6).   The  type  to  be  discussed  in 
this  paper  uses  the  free  precession  of  protons  as  described 
by  Bowen  (7)  and  by  Bacon  (8). 

The  free  precession  magnetometer  uses  a  strong, 

100  oersted,  magnetic  field  to  polarize  the  nuclear  moments 
thus  creating  a  large  value  of  magnetization  in  the  sample 

volume.   This  strong  field  is  suddenly  turned  off.  t   At 
this  time  the  nuclear  magnetization  vector  is  forced  by 
the  torque  of  the  earth1 s  field  and  the  intrinsic  spin 
angular  momentum  of  the  nucleli  to  precess  about  the  earth's 
field  vector.   The  frequency  of  precession  is  directly 
proportional  to  the  magnitude  of  the  earth's  field.  These 
precessing  nuclear  moments  induce  an  exponentially  decaying 
signal  in  a  detection  coil  surrounding  the  sample.   This 
signal  is  amplified  and  its  frequency  determined  in  order 
to  establish  the  value  of  the  earth's  field. 
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The  free  precession  system  has  a  limitation  in  the 
rate  at  which  measurements  can  be  made.   The  sample  must 
be  polarized  with  a  steady  magnetic  field  of  about  100 
oersteds  for  a  time  on  the  order  of  one  second.  Following 
this  there  is  a  frequency  count  time  which  consumes  another 
second.  With  careful  design  the  entire  cycle  can  be  com- 
pleted in  only  one  second  but  with  a  sacrifice  in  accuracy. 
Thus  the  data,  rate  is  limited  to  once  per  second.   For 
aerial  prospecting  where  the  velocity  over  the  terrain  is 
easily  100  yards  per  second  a  more  continuous  flow  of 
information  is  desirable  for  detailed  mapping.  Another  use 
for  a  high  information  rate  would  be  upper  atmosphere 
measurements  by  rockets  where  the  entire  time  of  flight 
for  several  hundred  miles  may  be  only  a  few  minutes. 

One  answer  to  obtaining  a  higher  data  rate  is  to  use 
two  systems  and  time  share  their  outputs  in  the  recording 
or  display  device.   To  do  this  two  major  things  must  be 
accomplished.  First,  two  systems  must  be  married  by  a 
proper  timing  cycle  and  the  dual  information  presented  in 
suitable  form.   Second,  the  polarizing  and  signal  coil 
of  each  system  must  be  magnetically  isolated  from  that 
of  the  other.   Each  system  has  one  coil  which  serves  both 
to  polarize  and  to  detect  the  signal.   The  isolation  is 
essential  since  the  polarizing  field  of  100  oersteds  at 
one  coil  must  not  affect  the  earth's  field  of  one-half 
oersted  at  the  other  coil.   This  is  so  because  one  coil  will 
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be  detecting  a  signal  at  the  same  time  the  other  is  being 
polarized,  A  distance  of  four  feet  was  selected  as  the 
coil  separation  distance  with  a  view  that  the  system  should 
eventually  be  airborne  and  probably  towed  in  a  "bird".   This 
paper  is  primarily  concerned  with  the  coil  designs  and  methods 
investigated  for  isolating  them. 

Suoh  a  dual  system  was  successfully  demonstrated  using 
a  double  concentric  solenoid  type  coil  for  each  system.   The 
timing  cycle  was  obtained  from  the  binary  counters  in  the 
frequency  determining  elements  of  a  standard  magnetometer. 
The  system  is  described  in  some  detail  in  Chapter  VI  and 
more  completely  by  Yowell  (14). 

The  units  of  measurement  employed  for  magnetio  fields 

are  often  puzzling  even  to  those  working  currently  in  the 
field.   Therefore,  the  following  comments  seem  in  order.   The 
geophysioists  and  most  workers  in  magnetism  have  found  it 
convenient  to  use  the  electromagnetic  units  (emu)  the  chief 
value  of  which  is  that  the  permeability  of  free  space  is  one 
gauss/oersted.   Therefore  numerically  the  magnetic  induction 
B  equals  the  magnetic  field  intensity  H  in  free  space  and 
for  practical  purposes  in  air.   At  the  International  Conference 
on  Physios,  London  meeting  of  1934,  it  was  decided  to  call  H 
by  the  name  oersted  and  B  by  the  name  gauss.   Gauss  was  formerly 
used  for  the  field  strength  H.   The  acceptance  of  "oersted" 


has  been  somewhat   less    than   overwhelming.      Consequently,   gauss 
is   still  used  frequently   in   the   literature  for   H.     Geophysi- 
oists  have   found   it  convenient   to  define   a  unit   of  magnetic 
field  strength  called  the   gamma  which  is    10"5   oersteds.      The 
earth's   field    is   thus  nominally  one-half  oersted   or  50,000 
gamma.      The   gamma    is    convenient  because   the    diurnal  variation 
of  the  earth's   field   is   on  the   order   of  a  few  tens   of  gammas. 

The   emu  system  is   based   on   the   idea  of  a  unit   magnetic 
pole.     H  is   measured   in   oersteds  where   each  oersted  represents 
a  dyne   of  force  acting  on  a  unit  pole   placed  at  the  point   of 
field   measurement.      In   the   meter  kilogram  second  system  H  is 
measured   in  units   called  ampere   turns  per  meter.      This   points 
out   the   fundamental  observed  fact  that  a  loop   of  steady 
current   exhibits   the   same   properties   as   the   pole   of  a 
permanent  magnet.      In  the  mks   system  the   permeability  of 
free   space    is   4  IT  x  10"7   henrys  per   meter. 

In  what  follows  both  mks  and  emu  units  will  be  used  in 
equations   and  suoh  will   be  noted  in  the   margin   if  ambiguity 
is   possible.      The   following   conversion  table   should  be 
helpful   in  accounting  for   the  4  7T's   and  powers   of  ten  which 
are  always   cropping  up   in  magnetio  problems. 


TABLE  I 

For  measurements  of  the  same  physical  quantity  the 
value  in  either  the  mks  or  emu  system  would  be  as  indicated 
below: 


SAME 

SYMBOL 

MKS  UNIT 

EMU  UNIT 

LENGTH 

I 

1 

METER 

102    CENTIMETER 

MASS 

M 

1 

KILOGRAM 

103  GRAMS 

TIME 

t 

1 

SECOND 

1        SEC  OND 

FORCE 

F 

1 

NEWTON 

105  DYNES 

WORK 

W 

1 

JOULE 

107   ERGS 

POWER 

P 

1 

WATT 

107   ERG/ SECOND 

ELECTRIC    CHARGE 

a 

1 

COULOMB 

lO"1  ABCOULOMB 

EMF,    POTENTIAL      # 

V 

1 

VOLT 

108   ABVOLTS 

CURRENT 

I 

1 

AMPERE 

10" 1   ABAMPERES 

RESISTANCE 

R 

1 

OHM 

9 
10      ABOHMS 

CAPACITANCE 

C 

1 

FARAD 

10~9    ABFARADS 

MAGNETIC    CHARGE 

m 

1 

WEBER 

108'47fUNIT  POLE 

MAGNETIC   DIPOLE- 
MOMENT 

m 

1 

WEBER- METER 

1010Mjr  POLE-CM 

MAGNETIZATION 

M 

1 

WEBER/ S^, 

■  MTR 

104'47T  POLE/SQ.CM 

MAGNETIC    FIELD- 
INTENSITY 

H 

1 

aMP. TURNS 

i/MTR 

4TTX10~3   OERSTED 

PERME  ABILITY 

/> 

47TX10"7    HENRYS /MTR 

1    GAUSS /OERSTED 

INDUCTION 

B 

1 

WEBER /SQ. 

,MTR 

104   GAUSS 

INDUCTION  FLUX 

? 

1 

WEBER 

108    MAXWELLS 

INDUCTANCE 

L 

1 

HENRY 

109   ABHENRYS 

CHAPTER  II 
NUCLEAR  MAGNETIC  RESONANCE 

1.  Theory 

Nuclear   magnetic  resonance    (NMR)   measurements  may  be 
possible    for  nuoleii  which  possess  a  finite   magnetic  moment^ 
and  an  angular  momentum  a.    a.  is  usually  expressed  as  a  spin 
I     according  to    the   relation        a-   =  I  h/2TT  2-1   emu 

where  h   is   Planck's   constant  «   6.6234  x  10~27   erg-sec. 
I  may  have   integral   or   half   integral  values.      The  value   of 
the  proton  spin   is   l/2.      For  raagnetio  field  measurements 

the   proton   is   usually  ohosen  as    the  nuclear  particle  to 
use.      It  has  a   comparatively  large   gyromagnetic  ratio 

fy    ~     ±ly"fl/lap\  2-2   emu 

( ju    and   a.    are   either  parallel   or   anti-parallel)   which 
makes  possible  a  highly  accurate   measurement   of   the   pre- 
cession frequency.     The  precession  frequency  is    given  by 
the  relation  CO  n  tfp  H  2-3  emu 

Electron  paramagnetic  resonance  has  been  used  for   earth's 
field  measurements  but  its   accuracy   is    limited  by   its 
broad  line   width  of  approximately  2000  gamma    (5). 

The  national   Bureau  of   Standards  has   measured    Yp  to   an 
aoouracy  of   one   part    in  40,000  as    (2.67528  ±    .00006)   x  104 
seconds"*1     oersteds"1.      Thus,    since    Yp     Is   known,   a  field 
can  be   accurately  and  absolutely  measured  by  making  a 
frequenoy  measurement  and  using   the   relation  of   equation 
2-3:      H  =        2TTf/*p    =   2.3486   x  10-4f(ops)  2-4   emu 
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To  exhibit  proton  resonance  it  is  necessary  to  have  a 
supply  of  hydrogen  to  provide  the  protons.   Water  is  a 
convenient  source.  The  sample  is  placed  in  a  large  magnetic 
field  HQ,  say,  10,000  oersteds.  As  a  simple  model  it 

can  be  assumed  that  each  proton  will  align  its  magnetic 
moment  axis  either  with  or  against  the  field.   There  is  a 

finite  energy  difference  W  between  the  two  positions  and  this 
energy  is  absorbed  or  given  up  during  a  transition  from  align- 
ment with  the  field  to  against  the  field  or  vice  versa.   This 
energy  is  equal  to  the  torque  exerted  by  the  magnetic  field 
on  the  magnetic  moment  of  the  proton  times  the  angle  through 
which  it  is  rotated,  namely  TK  It  is  apparent  from 

figure  2-1  that  the  torque  is  not  constant  but  varies  as 
sin  S    where  Q    is  the  angle  between  the  field  H0  and  the 
axis  of  the  proton.   Consequently: 

W   =.  ±  f  \\0yjp  £/„  Ode        =■  £  l/ip  Wo        2-5  emu 

From  quantum  theory  it  is  known  that  energy  is  absorbed 
or  radiated  in  discrete  amounts  equal  to  hf  where  f  is  the 
frequency  of  the  radiant  energy  and  h  is  Planck fs  constant. 
Thus  W  must  equal  hf  =  ZlSLyUp   or  H  =  ht/Z/fp   sl^f^:^,  2-3  emu 
Generalized  to  H  =  Ihi/u     this  equation  describes  NMR  for  all 
particles  when  I  and  u    are  those  of  the  particle  Under 
cons  ideration. 

As  another  approach  to  NMR  consider  a  sample  of  water 
which  is  in  no  external  field.   If  an  external  field  H 
is  suddenly  applied  consider  the  motion  of  one  of  the  randomly 


distributed  protons   whioh  make   an  angle  Q  with  H$   •      Since 
the  proton  possesses  angular  momentum    a   and  a  magnetic 
moment/W  ;    H0     will   exert  a  torque   on   the   proton  as   shown  in 
figure  2-2.      Torque   iB   equal   to   the   rate   of  change    of  angular 
momentum.      Therefore: 

£^    —/UfHo  Si*&  The  vector  <£  then  processes  about 

Hq     with  a  frequency,    say,   f.      The  linear  rate   at  which  the 
tip  of  a  moves  around   the   circle     of  radius    cx£>*£   is   ^    . 
Thus  |^  =2Tr{c*S.^^/0|+oS^a  Consequently  HQ  ■    ^£2>    =  °°/x>e 
which  is    equation  2-3  again.      It  must  be   pointed  out  that 
the  resonance   frequency   is    independent   of   the    angle    ©    . 

2.     Experimental   methods 

The  techniques  for   detecting  the  WMR  signal  are 
historically  two.      The  Harvard  group  headed  by  Furcell  used 
a  bridge   comparison  technique   to   detect   the    ohange    in  Q,  of 
a  Gravity  containing  protons  when  the  resonant  frequency  was 
applied    (10)    (11).      The   Stanford  method   of  nuclear   induction, 
directed  by  Bloch,   used  the  protons   to   couple   energy  from 
an  applied  radio   frequency  field   coil   to  a  pick  up  coil 
which  were   otherwise    essentially  decoupled    (12). 

In  the  nuclear   induction  method  a   10,000   oersted 
field  H0  is  applied  to  the   sample    in  the   Z  direction  and 
then  swept  a  few  oersteds  by  a   60  cycle   line   voltage    (13). 
As   shown   in  figure   2-3  a  radio  frequency  field  2H_oos6o  t 
at  42.5  megacycles   is   applied  by  small   coils    in  the   X 
direotion.     A  pick  up   coil  solenoidal  around  the   sample 
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Figure  2-1 
Proton  Alignment  i#  Steady  Field 


Figure  2-2 
Proton  Precession  in  Suddenly  applied  Field 
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with  axis  in  the  Y  direction  is  connected  to  a  narrow 
band  preamplifier  followed  by  a  homodyne  receiver.   The 
output  of  this  receiver  is  placed  on  the  Y  input  of  an 
oscilloscope  and  the  60  cycle  modulation  on  the  X  input. 
The  local  oscillator  signal  for  the  homodyne  receiver  is 
supplied  by  a  leakage  signal  from  the  transmitter  coil  to 
the  receiver  coil.   The  phase  of  the  leakage  is  adjusted  by 
means  of  copper  and  carbon  paddles  which  steer  the  leakage 
flux  through  the  pick  up  coil. 

As  the  H  field  is  swept  through  the  precise  value  of 
resonance  it  is  possible  to  have  displayed  one  of  two  types 
of  signals.   These  signals  correspond  to  the  proton  induced 
voltage  which  is  either  in  phase  or  in  quadrature  with  the 
applied  radio  frequency  signal.  Bloch  (12)  describes  these 
as  the  "v"  and  "u"  modes  corresponding  to  the  magnetization 
components  in  a  rotating  coordinate  system.  Figure  2-4 
shows  a  "v"  mode  signal  and  figure  2-5  a  "u"  mode  signal. 
The  sensitivity  associated  with  detecting  either  signal 
is  the  same.   Usually  the  symetric  or  TTv"  mode  is  displayed 
in  NMR  spectroscopy  work.  Figure  2-6  shows  a  portion  of 
the  ethyl  alcohol  spectrum.   The  field  difference  between 
the  major  peaks  is  a  few  milli-oersteds. 

3.   Signal  Voltage 

The  NMR  signal  is  not  induced  by  all  the  protons  in 
the  sample  but  rather  only  from  the  net  difference  between 
those  turning  to  the  lower  energy  state  and  those  turning 
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'V  Mode  Signal 
Figure  2-4 


"u"  Mode  Signal 
Figure  2-5 


Portion  of  Ethyl  Aloohol  Spectrum 
Figure  2-6 
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to  the   higher  energy  state.      This   fractional   difference    is 
given  by  the  Maxwell-Boltzmann  law  for  the   particle   distribu- 
tion according   to    energy  levels.      Let  n*      and  nh  represent 
the  number  of  low  and  high  proton  energy  levels  and  2/^H 
be   the  difference   in  energy  as   given  by  equation  2-5.      Then 
n£   «  nn     exp    (2/^H/KT)  2-6  emu 

Since   there   are   6.686  x  102^  protons    in  one   cubic   centimeter 
of  water  n  ■   n£  +  nh  *   6.686  x  1022.      Solve   these   two  equations 
for  the  difference  n^   -n^.     Assume    the   following  values 

for   the   parameters: 

yuf  -LH-IKIO'         ergs/oersted  Ho~  /0>oo°        oersteds 

/C  -/,32x/o        ergs/degree  Kelvin    '  -  ^00  Kelvin 

Therefore:      ^f>  H*/KT    =  &-Wo"^  <<  I 
Therefore:      e^P(^V^//<T)    -  I  -h  <Zyty  Ho//z  T 

and         n*-~"h   -   ^^//cf    =  Zi.'+Z  +  io'^     protons     2-7   emu 

The   ratio. of.    n£.-n*    shows  that   out   of   every  million  protons 

rv 

in  the  sample  only  three  contribute  to  the  net  signal 
voltage.   Fortunately  this  is  sufficient  for  detection.   In 
fact  much  smaller  samples  than  one  cubic  centimeter  and 
lower  fields  are  feasible. 

The  application  of  H  causes  an  excess  of  protons  in 
the  lower  energy  state  which  produces  a  macroscopic  magnet- 
ization  MQ  •  ;c  HQ  where  X  -    ' 3l»-  KT"   gauss/oersted  2-8  emu 
is  the  nuclear  paramagnetic  susceptibility.   For  protons 
in  water  A  =   3.4  x  lO"1^  gauss/oersted.   MQ  is  the  magnitude 
of  a  vector  quantity  which  has  i,  Y,  and  Z  components.   The 
signal  detected  in  the  NMR  apparatus  is  proportional  to  the 
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time  derivative  of  the  Y  component  of  the  magnetization. 

4*        ^t:  -5^-    abvolts     2-9  emu 

N  *  number  of  turns  in  the  ooil  and  A  is  the  sample 
oross  section  taken  perpendicular  to  the  coil  axis. 

M       —         Alp    Si*  Hit 

n*'~        U+S-y*  2-10 

2-11 


H       is    the  applied  steady  field  modulated  "by  the    60   cycle 
sweep.  Hr  »   tv/ty  Hrf  =   2^008  OJt:     defines    the 

applied  radio  frequency  field.      Bloch    (12)    shows   this 
solution  for   1^  is  valid  when  60    and  H       vary  adiabatically f 
that    is  when      ||~|  <<_|*,H,|  2-12 

At  resonance   H0  =   Hp,      E   =    0,    and  My  ■   MQ8inwt  £_j3 

This   results    in  V  «      4TT nAuj  M0  Coi  tot  k/d~*    volts  2-14 

Using   CO  ~fy  Hd         and    *vj  -    /£ff  gives: 


_    x  ^(^y/^^H*  a*u>t  A  /<T* 


V  "     ^FTcT  '^"*    ^  C<^     A/U  volts  2-16   emu 

Thus,   as    indicated  previously,    a  high  value    of   gyromagnetic 
ratio,   Y    ,    is   desirable   since   signal   voltage  varies   as     X~    . 
Except   for  tritium  the  proton  has    the   highest     o    known. 
The  value  for  tritium  is    only  seven  percent  higher  than   that 
for   the  proton.      Consequently  the   proton's   abundance   and 
high  gyromagnetic   ratio  make   it   a  natural   choice   for   field 
strength  measurements. 
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4.  Signal   to  noise  ratio 

The   thermal   noise   voltage   from  the   resistance    in   the 
detection   ooil   is: 

Taking  the  ratio   of    the   peak  voltages   gives: 

cy    -    LCZ^f/VA^2  Ho*  *\f*  2_18 

5.  Relaxation  times 

Another  .aspect   of   the   NMR  picture   is  the   matter  of 
relaxation  times.      Before  the   application  of  a  magnetic 
field  the  populations  are   equal    in   the    two   levels  but  the 
application  of   HQ  must  result    in  a  net  unbalance   in  favor 
of  the  lower  energy  level.      This   unbalance   builds  up  until 
an  equilibrium  is   reached  between  the  rate   of  absorption 
of  energy  by   the   protons  and  the   rate  at  which  they  are   able 
to   transfer   their   energy  to   the   thermal  reservoir.      The 
reservoir   is   comprised  of  the   other  motional   degrees   of 
freedom  of   the  particles   comprising   the    sample.      The   expon- 
ential  rate   at  which  this  equilibrium  condition   is   reached 
is   characterized  by  a  time  constant   T..  .      The   thermal 
reservoir   is   often  referred  to  as   the   lattice   and  T-,    as 
the   spin-lattice    or   longitudinal  relaxation  time.      T-^   may 
vary  from  a   few  hundred  micro-seconds   to  many  seconds 
depending   on   the  nucleous   and   its   environment. 

There    is   also  another   relaxation  time    T2   the   spin-spin 
or   transverse   relaxation  time.      This    is    caused  by   the   effect 
of   local   fields   at  the    proton  existing  beoause    of   the 
neighboring  protons.      The   same   observed   effect   is    also 
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caused  by   inhomogenieties   of   the   H     field.      This    is   apparent 
from  equation  2-3:   CO^Yft  .      If  H  is   not   everywhere    the   same 
in   the   sample    there  will   be   a  frequency  distribution  about 
6^0    s    *X    H0   .      The   time   for   these   individual   precession 
vectors   to  spread  out    so   the   signal  has   dropped  to  37$ 
is   Tg.      For  protons    in  water  T<?   is  also   about   two   seconds 
if  the   external  field  is    uniform. 
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CHAPTER  III 
FREE  PRECESSION  MAGNETOMETER 

1#   Need  for  free  precession  magnetometer 

The  free  preoession  magnetometer  is  necessary  for 

measuring  low  fields  because  an  NMR  technique  as  described 

above  is  not  practical.   The  signal  to  noise  ratio  available 

is  the  primary  limitation. 

Using  equation  2-18  and  the  typical  detection  coil 

parameter  values  given  below  results  in: 

S/N  s  10  H02 

N  =     500  turns  R  =  2     ohms 

A  «     100   square   centimeters  B  «   50  cycles 

Since   H0   is  the   earth*s   field  and  is    only  one-half  oersted 

the  resulting   signal   to   noise  ratio   is   2.5.      According  to 

Yowell    (14)   for  a  signal  to  noise  ratio  of  ten  to   one   the 

rms   error  for  a  one   second  oount  time  would  be   one-tenth 
gamma.      For  a   one- tenth  second  oount   time    it  would   be   five 
gamma.      In  order   that  a  system  have  an  accuracy  comparable 
with  other  magnetometers   the   error  should  be   about   one-tenth 
gamma.      It   is   apparent   then   that  a  one-tenth  gamma   accuracy 
and  a  one-tenth  seoond  data  rate   could  not  be   achieved  using 
the   NMR  technique. 

2.      Free   precession  magnetometer 

Russell   Varian    (6)    oonoeived   of   a   fundamentally   simple 
means   of    increasing   the    signal   several  hundred   fold  by 
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using  a  large,   100   oersted,   polarizing  field  whioh   then 
created  a  magnetization  roughly  200   times   as  large   as   the 
earth's   field  would   do.      In  equation  2-18  H  2  can  be  replaced 
by  HpHe  representing  the   polarizing  and  the   earth's   fields 
respectively.      Consequently  the   signal    to  noise   ratio  becomes 
500   to   one. 

When  this  D.C.   polarizing  field   is   applied   the   magnet- 
ization  increases   in  an   exponential   fashion  with   a  time 
constant  T,.      The  field    is   usually  left   on  for   one    or   more 
time    constants   and  then  cut   off  by  a   critically  damping 
RC  circuit.      When  the  polarizing  field  has   decayed  to  about 
ten  oersteds,    roughly  one-half  ampere,    the   circuit    is 
opened  completely.      The   fast  final   turn  off   is^essential    or 
else   the   proton  spins    oan  follow   the    decaying  field   and  no 
observable   precession   occurs. 

At   turn  off  the   protons   begin  a  precession  about   the 
earth's  field  and  the   coil   terminals   are  now  connected   to 
a  capacitor  which   tunes   the   coil   to  the    frequency   of  the 
expected  signal.      For   the  two   kilocycle   signal   corresponding 
roughly  to   an   earth's   field   of  one-half   oersted  this    signal 
will  build  up   in   the  resonant  circuit  with   a  time    constant  = 
QJu  .      For   the  normal  Q  of   80   this    is   6.5   mill  i-seconds . 
The   signal    itself  will   decay   in  an  exponential    fashion  with 
the    characteristic   time   constant   Tp.      For  water  T^   is    about 
two  seconds.      Thus    if  the    initial   signal    to   noise   ratio  was 
500   to    one   and  a   count   is   made   for  two   time   constants    the 
signal   to  noise   ratio   at   the   end  of   the    count   will  be 
S/N  =   500C2   =    67 
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This    Is  well  within  limits   for   the  prescribed  accuracy   of 
one-tenth  gamma.      Figure  3-1  shows   a   free   precession  signal. 

It  has  been  assumed  that  the   earth's   field   is   always   at 
right  angles   to   the  polarizing  field  for   this   provides   the 
maximum  signal.      If  the   two    fields  were   aligned  parallel 
no   signal  would  result.      In  fact   the  signal   magnitude    is 
proportional  to  the   sine   of   the   angle  between   the   polarizing 
field  vector  and  the    earth's   field  veotor. 

The  free   precession  magnetometer  measures   the 
earth's   magnetic  field  by  measuring   the   proton  precession 
frequency.      To  do   this    there   are  basically  two   functions   to 
be   performed.      One    is    the  polarizing  of  the  sample   and   the 
other    is    the   frequency  determination.      The   frequency  deter- 
mination could  be   made   directly  with  a  counter  or   on  a  beat 
frequency  basis  with  a   standard  oscillator.      However,    in 
counting  there   is   always  a  possible   error  amounting  to  one 
count.      For   example,   a  discrete   counter  could  register 
either  2000  cyoles   or  2001  cycles   but   the   actual   frequency 
might  lie   anywhere  between  these   two  numbers.      Thus   for  a 
two   kilocycle   signal   there  would  be   a   one    in  2000  error  or 
an  error  of   25  gamma   in  an  earth's   field   of  50,000  gamma. 
To   reduce   this  error  the   two  kilocycle   signal    is   used  to 
gate   a  high  frequency  crystal   driven   counter.      If  the 
crystal  frequenoy   is  500KC   the    least  count   error   is   reduced 
by  a  factor  of  250   or   to   one-tenth  gamma. 

Figure   3-2   is   a  block  diagram  of   the   Varian  V-4900 
station  magnetometer.      Figure  3-3  shows   the   packaging   of 
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FIGURE  3-1  FREE  PRECESSION  SIGNAL 
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Station  Magnetometer  Blook  Diagram 
Figure  3-2 
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FIGURE      3-3 
STATION  MAGNETOMETER 


the  magnetometer  and  figure  3-4  shows  two  of  the  detection 
coils,  a  brief  description  of  the  station  magnetometer 
follows. 

The  coil  is  wound  with  five  layers  of  97  turns  each 
on  a  4.25  inch  outside  diameter  coil  form  six  inches  long. 
The  wire  is  #15  formvar  coated.   This  results  in  an 
inductance  of  15  millihenries  and  a  Q,  of  80  at  two  kilocycles 
There  is  a  75  foot  twinax  cable  connecting  the  coil  to  the 
polarizing  relay. 

The  polarizing  relay_ is  a  double  pole  double  throw  D.C. 
actuated  type.   Two  contacts  are  necessary  since  the  coil 
terminals  are  balanced  to  ground  while  one  side  of  the  polar- 
izing supply  is  grounded.   The  switch  has  an  RC  circuit 
across  it  providing  an  overdamped  current  decay.   When  the 
polarizing  current  has  dropped  to  about  one-tenth  of  its 
original  value  of  five  amperes  the  field  has  decreased  to 
about  ten  oersteds  and  a  relay  completely  opens  the  circuit. 
To  be  effective  the  final  turn  off  must  occur  in  a  time  much 
less  than  the  period  of  the  two  kilocycle  precession  signal. 
Consequently  the  self  resonant  frequency  of  the  coil  should 
be  20  kilocycles  and  preferably  higher. 

The  detection  coil  is  series  tuned  and  transformer 
coupled  to  the  first  of  four  stages  of  audio  amplification. 
Following  amplification  the  signal  is  squared  and  fed  to 
the  first  of  thirteen  flip-flops  comprising  the  signal 
counter.  When  enabled  by  a  signal  from  the  timing  cams 
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the  binaries  begin  counting  from  their  reset  value.   This 
reset  value  is  not  zero  but  effectively  a  negative  number 
corresponding  to  128  cycles  of  the  signal.   This  eliminates 
the  transient  rise  of  the  signal  amplitude  from  the  part 

which  is  counted.   When  the  binaries  reach  zero  the  500  KC 
gate  opens  beginning  the  count  of  the  500  KC  oscillator  signal. 
This  gate  is  closed  stopping  the  500  KC  count  when  the  pre- 
cession frequency  has  completed  256,  512,  1024,  2048  or 
4096  cycles.   The  count  interval  is  selected  by  switching 
out  successively  the  last  five  binary  stages  in  the  signal 
counter. 

The  crystal  oscillator  is  an  electron  coupled  Pierce 
oscillat.or  running  at  100  KC.   The  output  of  this  is  amplified 
and  multiplied  by  five  in  a  second  stage.   This  500  KC  signal 
actuates  a  squaring  amplifier  which  is  gated  to  the  high 
speed  binary  counters. 

There  are  nine  sets  of  binaries  in  the  500  KC  counter 
so  that  a  count  of  29  or  512  is  possible.  At  the  end  of  the 
count  the  bincry  number  is  analogued  as  a  DC  voltage  level 
by  varying  the  current  through  an  analoging  resistor.   The 
current  is  varied  by  introducing  any  of  nine  different 
current  paths  through  relays  selected  by  the  state  of  the 
binaries.   Only  the  last  512  cycles  or  fraction  thereof  in 
the  500  KC  count  is  recorded.   If  the  signal  count  is  made 
for  2048  cycles  of  a  2000  cycle  per  second  signal  then  this 
is  (2048/2000)  500,000  *  512,000  cycles.   The  record  then 
is  made  of  only  a  fraction  representing  the  last  one  thousandth 
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part  of  the  count.  Using  H  »  ZRf/ltp     to  determine  the  field 
corresponding  to  a  2000  oycle  per  second  signal  gives 
H  s  2  x  ,23486  -  .46792  oersteds.   The  recording  range 
represents  a  thousandth  part  of  this  or  46.8  gamma.   By  switch- 
ing out  the  last  binary  stage  and  doubling  the  analog  resistor 
the  recorder  sensitivity  is  increased  to  23.4  gamma  full  scale. 
The  most  sensitive  scale  of  11.7  gammas  is  obtained  by  dropping 
the  last  two  binaries  and  again  doubling  the  analog  resistor. 

The  recorder  is  a  standard  Gr-10  graphic  recorder  made  by 
Varian  Associates.   The  plot  is  rectangular,  full  scale  five 
inches,  with  a  chart  speed  varying  from  4  inches  per  minute 
to  four  inches  per  hour.   The  analog  voltage  is  only  avail- 
able from  the  end  of  the  count  until  the  time  of  reset; 
therefore  the  recorder  servo  is  supplied  with  power  only 
during  this  time. 

The  timing  cam  unit  consists  of  a  60  rpm  timing  motor 
and  five  disc  cams  actuating  micro-switches .   The  functions 
performed  are  actuating  the  polarize  relay,  triggering  the 
signal  count  gate,  resetting  both  counters  and  turning  on 
the  servo  motors  in  the  recorder. 

Figure  3-5  shows  a  portion  of  a  station  magnetometer 
record. 
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Figiire  3-5   Station  Magnetometer  Record 


CHAPTER  IV 
COIL  DESIGN  AND  CONSTRUCTION 

1.  General 

The  high   information  rate   magnetometer  requires  that 
the   earth's   field  be   measured  at   two  points   four   feet   apart 
using   essentially  two   station  magnetometers   and  time   sharing 
their   outputs    in   the   recorder.      This    means   that  to  have   an 
aoouracy  on  the   order   of   one-tenth  gamma   the   polarizing 
field   of  one    coil   should  not   contribute   more   than  one-tenth 
gamma    in  a  vector  addition   to   the    earth's   field  at   the   other 
coil.      At  first   the   problem  seems    easily  solved  by  aligning 
the   residual   field  from  the    coil   at   right   angles   to   the 
earth's   field.      However,    this  would   do  only  for  a   stationary 
system  or  a  moving  system  properly  stabilized.      To  have   to 
stabilize   the   system  would   sacrifice   the   most  advantageous 
feature   of   the  free   precession  system.      That    is   that  the 
measurement   is    independent   of  coil  angular  relation  to 
the    earth's   field  vector. 

The   problem  involved  here   has   nothing   to   do  with  the 
mutual   coupling  between  the   coils.      The  effect  which  must 
be   reduced   is   the  vector   addition   of    the   D  C  magnetic  polar- 
izing field  of   one    coil   to   the    earth's   field  at   the   other 
oo il.      However,    if   there    is  a   ripple    in   the    polarizing  power 
supply  then  the   mutual    inductance   will   cause   a  ripple   noise 
modulation   of   the   signal.      Consequently  low  mutual    is 
desirable    and   is   achieved  as   a   natural    result   of  solving 
the   prj^ary  problem. 
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2.  Design  Parameters  and  Restrictions 

Any  coil  designed  for  this  system  should  attempt  to 
obtain  a  high  signal  to  noise  ratio.   Equation  2-18  gives 
the  signal  to  noise  ratio  and  from  this  the  parameters 


which  influence  coil  design  may  be  obtained. 
% 


rv     __     A  //ArtK£+0*S"  *M$z  4-1 


Where  H  He  has   replaced  HQ2   as   explained  before. 
Eliminating  constants   and  factors  which   are   not   affected 
by  coil  design  results   in: 

Where:   n   (number  of  nuoleii)  oC    Volume       -JiA    <=<-  £d 

0 

JL    *  ooil  length  A  =  ooil  diameter 

B  (band  width)   «%   s      %^\^  <   % 

Therefore:      S/N     oc     NA     ^^  4-3 

Once   a  wire  size   is   chosen  R  becomes   practically  proportional 
to  K  and  if   JL    is  chosen  equal   to  d  then: 

S/N    oC    A/^d  4-4 

This  dearly  indicates  the  bigger  the  coil  the  better. 
Praotioally  there  are  limitations.   In  arriving  at  equation 
4-4  the  sample  diameter  and  coi^  diameter  were  assumed  the 
same.   Practically  this  can  not  be.   In  fact  using  large 
wire  to  reduce  the  resistance  means  there  will  be  a  consider- 
able volume  of  copper  which  can  not  be  used  for  sample. 
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A  better  approach  to  coil  design  would  be  optimizing 
the  various  parameters:  wire  size,  coil  length  and  diameter, 
and  number  of  turns  for  maximum  S/N  ratio.   This  to  be  done 
given  a  prescribed  volume,  or  weight  or  polarizing  power. 
In  airborne  or  rocket  work  all  three  of  these  restrictions 
would  be  important.   Since  optimum  coil  design  does  not 

bear  directly  on  producing  a  low  external  field  coil  it 
will  not  be  pursued  further  here. 

Because  the  polarizing  current  turn  off  must  be  rapid 
there  is  a  very  important  restriction  on  the  coil  design 
requiring  that  it  have  low  distributed  capacity.  At  least 
it  must  be  small  enough  so  that  its  self  resonant  frequency 
is  20  KC  or  greater.   The  exponential  decay  of  the  current 
should  be  rapid  -  at  least  a  fractional  millisecond  time 
constant.   The  decay  is  more  rapid  the  higher  the  resistance 
of  the  coil.   Consequently  the  signal  to  noise  and  the  turn 
off  requirements  are  opposite  in  their  demands  on  coil 
resistance. 

Increasing  the  Q,  of  the  coil  will  increase  the  signal 
to  noise  ratio.   However,  there  is  a  practical  limit  imposed 
by  the  diurnal  field  fluctuations  for  a  stationary  system 
or  by  the  field  contours  traversed  by  a  moving  system.   If 
the  Q,  is  100  then  the  pass  band  at  two  KC  is  20  cycles 
corresponding  to  500  gamma.   A  nominal  value  for  the  North- 
South  gradient  of  the  earth's  field  is  five  gammas  per  mile. 
Consequently  an  aircraft  could  easily  fly  beyond  the 
tuning  range  provided  by  20  cycles.   For  aircraft  systems 
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automatic}  frequency  tuning  by  bands  is  essential.   The 
dual  coil  system  discussed  herein  did  not  incorporate 
such  tuning. 

In  order  to  reduce  the  field  interaction  between  the 
two  coils  there  are  three  adjustable  variables.   First, the 
separation  distance  which  has  arbitrarily  been  taken  as 
four  feet.   Second,  the  polarizing  current  can  be  reduced 
thereby  reducing  the  field.   This  is  possible  if  the  signal 
to  noise  ratio  can  be  maintained  high  enough  by  other  means. 

Third,  by  configuring  each  coil  in  a  particular  manner  the 

f 
external  field  can  be  minimized.   Coil  configuration  is  the 

primary  technique  considered  in  the  remainder  of  this 

chapter. 

3.      Magnetic   Field   Calculations 

In  order   to  appreciate   the  magnitude   of  the  values 
in   the   coil  design  problem  it   is   necessary   to  digress   and 
solve   two  magnetic  field  problems.     First,    the    field   from  a 
wire   of  finite   length  carrying  a   steady  current   is   given  by: 
H  =    i^rA-2  f     ^£m9ik  From  figure  4-1  4-1   mks 

9-7T  J&t  ^Jl.  & 

Therefore        H  »  -^=~z    j    Si*Q<k9  amp- turns /meter  4-2  mks 

Solving  this   for  a  one   inch  length   of  wire  with  a,  equal   to 

four   feet  gives:        H  =    !?£_  Coi€^3i     -    7JI7lrt('0{  °^ 

_j 

since     Cbtg.§  -  2*9-8  -  .0/02    -  C»s9,  Therefore: 

TT  -3  / 

H  «  1.33  I  x  10    amp. -turns /meter  «  1.67  I  gamma 
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Figure  4-1 
Magnetic  field  from  a  Straight  Wire 


H-ft 
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Figure  4-2 


Magnetic  field  from  two  Straight  Wires 
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Thus,   a   current  of  one   ampere   in  a  single   one   inch  lead 
oan  cause   a  large   field   in  relation  to    the   one-tenth  gamma 
field  desired. 

The  second  problem  is    to   see  what   field  results   from 
two   such  wires   placed   side  by  side  but  carrying  current   in 
opposite   directions.      Taking   #14  wire  which  is    .0641   inches 
in  diameter  and  using   superposition,   a   figure   eight  pattern 
results  which   has  a  maximum  in  the  plane   containing  the 
wires.      This   maximum  is   given  by: 


H  *  JL     cos©      (    1'22 -        1*22 \      1.341x10        fn.  & 

zrr     oob^     ^1#22  .    .000815         T7ZZ .00081b) ~  ~Trf ' 

In   this  case   the    interest    is    in  the   effect   of   the    lead 

in  cable    to  the    coil.      If   its   effective   length   is   four 

feet   then  tan  ©,    =   2  and  cos  0,    *      .446  and 

-r 

H  mTrr    x  ,446  x  1#34  x  10"3  =  9.51  x  10"5  amp- turns /meter 
§r  H  s  .119  I  gammas.   Consequently  it  is  desirable  to 
make  the  leads  coaxial  as  this  will  reduce  the  external 
field  to  zero. 


4.  Coil  Configurations 
A.  Double  Solenoid 

There  are  several  possible  configurations  to 
achieve  a  coil  system  which  has  a  low  external  field.  As 
configuration  that  would  achieve  this  is  two  solenoids, 
one  inside  the  other.   If  the  flux  from  each  coil  is  the 
same  in  magnitude  but  opposite  in  sign  then  the  external 
field  will  be  small.   The  residual  field  is  due  to  the 
fact  that  the  shape  of  the  field  is  not  quite  the  same 
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sinoe   the   coil   sizes  are   different.      The   field  at   the   center 
of  a   long   solenoid   is   given  by     H  =     NI/1  and  the 

total   flux  by   <f  -  5 A   ~/*HA  -  <*"IA/i 
Where   N  is   the    total   turns   and  £     is   coil   length.      For  the 
total   flux  to  be   the  same:        N     I     A*/g        =     NglgAg/tf  p 

For   praotical  reasons    the    current  must  be   the  same    in  both 
ooils   and  X-^   is   chosen  equal   to    ^£    •      This  requires  N^A^«   NgAg 

or    /V/  <*,     ~Mzo<l     where    Ol   is  the    coil  radius. 
The  field  strength  at  the    center  of  the    concentric   coils 
is  H  »    -j—     (N^  -Ng)      ,      Therefore   there    is   a  net  field   inside 
while  the   external   field   is   greatly  reduced. 
B.   Flux  Ball 

There   is   a  configuration  for  which  the   external 
field  can  theoretically  be  made  zero.     This    is   the  double 
"flux  ball"  arrangement.      If  a  spherical   shell  has  a 
surface  current  distribution  such   that  the   flow   is   along 
parallels   of  latitude  and  varies  as  sine  Q   where  Q  is 
the   oolatitude  angle    then  the    field   inside    is   uniform. 
Further,   according   to  Blewett    (15)    the    external   fields 
are   described  in  cylindrical   coordinates  by: 

U        -    Q 

n^  a  *  sphere  radius 

Provided:   I(f  =  ^£_  =  **£  S>^& 

or  H,   =  ^^-^- 
Consequently  if  Ip  varies  as  %^   then  H   is  constant. 
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If  two    "flux  balls"  are   constructed  and  made   concentric 
then  the   coordinates    (^;<p;£)   are    the  same   for  both. 
Consequently    H^>,     can  be   made   equal  to  H^j.        if  Ho^z\\0  ax      4-3 
A  similiar  requirement   obtains  for   H_.     Practically  in 
order   to  vary   I y    the   number  of   turns    is  varied  along  the 
surface   of   the   sphere   in  accord  with   sine  0     .      If  N   is   here 
used  for  the   turns  per  centimeter  at  the   equator   then  to 
satisfy  equation  4-3   it   is   necessary  that    A/,  <*,     ~h/zCKx 
This   configuration  does   supply  two   external   fields  whose 
shapes   are  the   same   and  perfect  cancellation  would  be 
possible    if  perfect   construction  were  feasible.      Time   and 
expected  construction  difficulties   precluded  making  this 
type  coil. 

C.     Toroidal   Coil 

A  long  solenoid   bent   into  a  circle   makes   a   toroid. 
Under  the   assumption   that   the  bending  does  not   disturb 
the  solenoidal   windings   the  flux  through   the    toroid  should 
all   remain   internal    to    the   winding.      In  practice   this 
perfection   is   not   achieved.      However,    if   it   could  be   closely 
approached   then  a  low   enough  external   field  for   the   purpose 
in  mind  might   be   obtained,      a  toroidal   coil  has   a   distinct 
advantage   over  the   double    solenoid   or   the   flux  ball 
arrangements    in   that  no   matter  how  it*  is   oriented  with 
respect   to   the   earth's   field  approximately  half  of   its 
sample   volume  will  be   contributing   to   the   signal.      This 
is  not  so  for   the   solenoid  or  flux  ball   since   no  signal 
is   induced   if   the    polarizing  field  and  the   earth's   field 
are   parallel. 
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5.  Double  Solenoid  Coil  Construction 
A.  Double  Solenoid  #1 

The  first  ooil  system  constructed  was  a  double 
concentric  solenoid.   The  coils  were  made  so  that  their 
magnetic  moments  were  equal  by  making  /V,  <a,   -  /Vt^-a. 
This  coil  system  was  made  with  a  standard  station  magnetometer 
ooil  as  the  inner  solenoid.   This  coil  has  five  layers  on 
a  4.25  inch  outside  diameter  ooil  form  totaling  485  turns. 
Since  the  number  of  layers  was  odd  the  leads  did  not  terminate 
at  the  same  end  and  it  was  necessary  to  add  a  sixth  layer 
of  eight  turns  in  a  high  pitch  helix  to  bring  the  leads  out 
the  same  end.   The  outer  coil  had  4  layers  of  89  turns  each. 
This  was  known  to  be  more  turns  than  necessary  on  the  outer 
ooil  but  experimental  determination  of  the  exact  number 
seemed  called  for.   The  formvar  was  scraped  from  a  strip 
down  the  outer  layer  and  a  sliding  wire  contact  made  so  that 
the  proper  number  of  turns  could  be  determined.   Thirty 
turns  on  the  outer  layer  proved  adequate  to  reduce  the  field 
to  about  ten  gamma  per  ampere  at  a  distance  of  four  feet. 
It  seemed  that  the  lead  in  wiring  and  lack  of  symetry  were 
causing  most  of  this  field.   The  outer  layer  was  taken  off 
and  rewound  as  a  30  turn  helix  of  even  pitch.   Some  improve- 
ment was  noted  here  but  not  an  appreciable  amount. 
B.   Double  Solenoid  #2 

In  view  of  the  relatively  large  effects  of  short 
leads  and  winding  errors  a  second  concentric  solenoid  pair 
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was  constructed  with  six  layers  on  a  three  inch  outside 
diameter  coil  form  and  two  layers  on  a  five  and  25/32  inch 
coil  form.  Each  layer  was  six  inches  long  and  had  88  turns 
of  number  14  formvar  coated  wire.  Figures  4-3  and  4-4  show 
the  construction  of  this  coil.   The  inductance  of  the  coil 
was  7.1  millihenries  and  the  self  resonant  frequency  was 
40  KC.   The  field  from  this  coil  was  .36  gammas  per  ampere 
at  a  distance  of  four  feet  at  a  minimum  point  in  the  field 
pattern.  Figures  5-3  to  5-5  show  plots  of  the  field  compon- 
ents at  points  about  the  coil. 
G.   Source  of  Signal 

The  field  plots  have  been  given  in  gammas  per 
ampere  because  the  polarizing  current  is  a  variable  that 
can  be  changed  to  reduce  the  external  field.   Changing  the 
current  will  change  the  polarizing  magnetization  and  thereby 
the  signal  to  noise  ratio.   The  integral  of  H  over  the 
sample  volume  is  the  best  measure  of  the  effective  magnet- 
ization created  by  any  coil  configuration.  For  the  double 
solenoid  it  is  apparent  that  there  are  two  volumes  which 
can  contribute  to  the  signal:  the  inner  volume  and  the 
volume  between  the  coils.   The  field  is  greatly  reduced 
external  to  the  two  coils  so  this  volume  is  neglected.   The 
central  volume  field  is  about  100  oersteds  as  the  result 
of  the  two  fields  bucking.  The  field  between  the  coils 
is  the  resultant  of  adding  the  outer  coil  field  to  the 
portion  of  the  inner  coil  field  that  threads  the  space 
between  the  coils.   This  field  value  is  about  75  oersteds. 
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It  is  logical  now  to  determine  if  the  two  sample  volumes 

contribute  signals  which  are  in  phase  and  thus  enhance 
the  signal. 

A  double  subscript  notation  is  convenient  to  represent 
the  various  voltages  induced.   The  first  subscript  will 
indicate  the  sample  volume  producing  the  signal  and  the 
second  subsoript  which  coil  it  is  induced  in.   Thus  VAg  is 
the  voltage  induced  in  the  outer  coil  (2)  by  the  protons 
in  the  inner  sample  (A).   The  protons  are  initially  pro- 
cessing in  planes  parallel  to  a  plane  containing  the  center 
line  of  the  coils.   Since  the  fields  in  the  two  volumes  are 
in  opposite  directions  and  the  coils  are  connected  in 
series  opposing  the  following  equation  for  the  signal 
voltage  results:  V  -  VA1  +    vB1  +     Y^2     -  VA2         4_4 

Using  both  volumes  results  in  a  net  signal  gain  over  the 
signal  from  the  inner  volume  alone. 

In  order  to  use  the  outer  annular  volume  it  was 
necessary  to  insulate  the  windings  from  the  water.   This 
proved  to  be  a  difficult  task.  A  commercial  resin  "Epon", 
was  used  to  cover  the  coil.   This  proved  unsatisfactory 
as  after  a  few  days  a  leakage  path  of  one  magohm  appeared 
between  coil  and  water.   Since  the  system  required  two 
double  solenoids  the  second  one  was  insulated  with  Epon 
and  a  liquid  latex.   This  proved  satisfactory  for  the  short 
time  tested.   However,  the  final  system  tests  were  made 
using  kerosene  as  a  sample  since  it  is  an  insulator. 
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6.      Toroidal   Coil   Construction 

To   test   the   toroid  idea  a  toroidal   coil   form  was   made 
by  joining   two   "return  U"   pieces   of  three    inch  diameter 
pyrex  glass   tubing   into  a   circle.      The   wire  was   then  rolled 
on  by  leaving  the  wire   end  free.      Each  layer  was  pulled 
down  tight   and  cemented  with  Epon  before   the   next   layer 
was  wound.      The   wire   ends  were  welded   together  by  the 
hell-arc  process.      The   layers  were  wound  back  and  forth 
rather   than   continuously  around   in  order  to   avoid  creating 
current   loops  around  the   torus. 

If  the   radius   of  a  section   of  the   torus   is   "a"  and 
the   distance   from  its   center   to   the   torus   axis   is    "b"   then 
the   following  formulas  apply: 
Volume     =        Zir%c£    -  3-06%  Liters 

H  center  s   -~-£  =     15.1   I     oerstedB 

I     •/*>  A/*B>-(bx-cftS}         ■     9.5  mh. 
fhAV   =  TTNo^I  s 46  5*00 1      oersted  -   cm3 
The  numerical  values   are   those   of  the  constructed  four 
layer  toroid  having  856  turns.      The   integral   of  the   H  field 
over   the  volume    is    comparable   to   that   of   the   station 
magnetometer   coil   nominal   value   of  40,000   I   oersted-cm   . 
The   least  field  at   four  feet  was   two   gammas   per  ampere. 
The   Initial   signal    to  noise   ratio  was   72   to   one   for  a   six 
ampere   polarizing  current   and  a  3.3   second  polarize   time. 
It   is    apparent   that   this   coil  will  not  produce   the  desired 
low  field  unless   a  much   smaller   ourrent    is   used. 

In   order   to    increase   the   signal   to   noise   ratio   two 
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more  layers  were  added  giving  a  total  of  1270  turns. 
Figure  4-5  shows  this  coil.   The  inductance  of  this  ooil 
was  21.4  millihenries.   The  least  field  was  three  gammas 
per  ampere  at  four  feet.   The  signal  to  noise  ratio  was 
180.   Thus  there  has  been  some  gain  since  the  external 
field  has  gone  up  only  50$  while  the  signal  to  noise 
ratio  has  increased  by  150$.   Consequently  less  polarizing 
current  can  be  used  for  the  same  signal  to  noise  ratio 
and  thereby  the  external  field  is  reduced. 

This  six  layer  toroid  had  one  major  drawbaok.   Its 
distributed  capacitance  was  high.   The  ring  frequency  was 
13  KC  and  the  decay  time  constant  was  long  -  one  milli- 
second.  The  effect  of  this  was  to  oause  an  erratic  polar- 
izing current  turn  off.   This  was  evidenced  by  a  large 

(25$)  fluctuation  in  the  initial  amplitude  of  the  signal. 

i 

The  toroid  did  not  have  as  low  an  external  field 
as  the  double  solenoid  so  it  was  not  used  in  the  system 
test.  A  toroid  wound  on  a  form  which  did  not  have  such  a 
large  torus  diameter  as  three  inches  would  have  less 
winding  error  and  more  closely  approximate  an  ideal  toroid. 
Sample  volume  would  be  sacrificed  but  for  a  given  ourrent 
the  external  field  should  be  less. 

7.   Lead  in  Cables 

To  minimize  the  lead  in  cable  field  a  double  coaxial 
cable  was  constructed  by  sliding  a  one-half  inch  braid  over 
a  standard  RG-58/U  coaxial  cable.   The  center  wire  and  the 
inner  coax  carry  the  polarizing  ourrent  and  the  outer  braid 
is  tied  to  the  Faraday  shield  around  the  coil. 
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In  dealing  with  cables    it   can   be   pointed  out   that 
some    cables   available   have   steel   strands   in   them  for   strength 
and  care   must  be   taken  not    to  use   these.      Even  this  small 
amount  of   magnetic   material  near  the   coil  will   cause   field 
inhoraogenieties   that   can  measurably  reduce   the    observed  Tp. 
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CHAPTER  V 
COIL  FIELD  MEASUREMENTS 

1,  Measurement  Methods 

The  magnetic  field  about  the  constructed  coils  had 
to  be   determined.      The  primary  method  used  was   to   pass   a 
400  cycle  current   through  the  magnetometer   coil  and  then 
measure  the   induced  voltage   in  a  small   10,000  turn  pick  up 
coil   placed   at  a  measurement   point,     "As  a   check  on  this 
measurement  a   standard  station  magnetometer  was  used. 
The   magnetometer  recorded  the   change    in  field   occur ing 
when  a   known   DC   current  was   passed   through   the  low  field 
coil  placed   in  proximity  to  the   station  magnetometer   coil. 
Figure   5-1   shows   an  example    of   one   of   these   measurements. 

2.  Pick  up  Coil   Calibration 

The  10,000  turn  pick  up  coil  was  wound  on  a  two  inch 
diameter   split  ring  brass   coil   form  as    shown   in  figure   5-2. 
The   depth  of  winding  was    .625   inches.      Using  a  formula 
suggested  by  Frowe    (16) 

a(eff)s    [ao+a^Wj^ 

the   effective  radius   of   the  coil  was    .0337  meters. 

The  voltage   induced   is   then: 

Vs   Hjf  =A/FcL-d&    -/VircL^cuHcc^Lot  5-1   mks 

Where  H  sinujt     is   the  uniform  field  threading  the    coil 

at  right  angles   to    the  plane   of  the    coil. 

V  *   Nrra^aTr-fH    =  )o*TT(.o337)\7r»to7»2nf  H  -  a.32f  H«i6"^  5.2   mks 

Where   V  is    in  volts   and   H  is   in   ampere    turns   per  meter. 
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Converting  to   oersteds  results   in: 

V  volts   =  3.«  »  $  H0*iM    =  -2^^^VHflmteis  5-3 

H  oersteds  s  iHL-L  * Yvoits  5-4 

T 

H  gamma  -       t^9.    y.   Vmillivolts  5.5 

To  verify  this  calibration  a  pair  of  Helmholtz  coils 
of  17  inch  radius  and  17  inch  separation  were  used  to 
experimentally  calibrate  the  pick  up  coil.       —  ■*» 
McComb  (1),  page  132,  gives  the  following  convenient 
formula  for  the  field  at  the  center  of  a  Helmholtz  pair. 

H  =   89.9  N  I  5_6 

Where  H  is  in  gamma,  I  is  in  milliamperes ,  "a"  is  the 

coil  radius  and  separation  in  centimeters  and  N  is  the 

number  of  turns  on  each  coil. 

The  pick  up  coil  was  placed  at  the  center  of  and 

parallel  to  the  Helmholtz  pair.  A  400  cycle  current  of  a 

few  milliamperes  was  sent  through  the  Helmholtz  coils  and 

the  output  voltage  of  the  pick  up  coil  measured.   The  mean 

of  several  measurements  gave  3.27  millivolts  induced  in  the 

pick  up  coil  for  each  milliampere  of  current  in  the  Helmholtz 

coils.  Using  equation  5-6  and  the  values  of  the  particular 

Helmholtz  pair  results  in: 

H  gpmma  *  35.4  IM  »  J>^     W    *  10.84  W       5-7 

This   is  at  400  cycles.     The  voltage    induced   is   directly 

proportional    to   the   frequency  consequently: 

H  gamma  ■      10.84  x  400        Vmv  «  4340        mv_  5.8 

f  f- 
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This  compares  within  3$  of  the  value  given  by  equation  5-5. 
Experimental  error  in  this  calibration  was  approximately  2$. 

3.   Field  Measurements 

The  field  plots  of  the  various  coils  as  shown  in 
figures  5-3  to  5-8  were  made  using  the  10,000  turn  pick  up 
coil  connected  by  a  coaxial  oable  to  a  Hewlett  Packard  400D 
voltmeter.  With  this  arrangement  the  normal  noise  voltage 
was  ,04  millivolts  rms  corresponding  to  .447  rms  gamma  by 
equation  5-8.   This  indicates  that  a  one  gamma  field  could 
be  measured  and  the  error  due  to  noise  would  be  10$. 

There  are  two  ways  to  increase  the  signal  to  noise 
ratio  so  that  measurements  of  lower  fields  are  possible. 
The  first  is  to  increase  the  current  in  the  coil.   An  110 
volt  400  cycle  power  supply  was  used.   For  the  double 
solenoid  this  gave  a  current  of  6.2  rms  amperes  and  2.06 
rms  amperes  for  the  six  layer  toroid.  For  most  of  the 
measurements  out  to  four  feet  this  was  adequate.   A  second 
method  to  increase  sensitivity  is  to  take  the  output  of 
the  voltmeter  and  feed  it  to  a  400  cycle  band  pass  filter 
and  present  this  on  an  oscilloscope.  By  this  means  a  one- 
tenth  gamma  measurement  was  possible.   Since  the  radial 
fall  off  of  the  field  was  in  general  of  the  form  H  =  Kr"x 
it  was  possible  to  make  measurements  at  three  feet  and 
calculate  the  value  at  four  feet. 

The  field  components  were  measured  in  cylindrical 
coordinate  fashion  by  placing  the  pick  up  coil  so  the  -o; 
f »  or  Z  coordinate  vector  would  thread  the  coil.   The 
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External  Magnetic  Field  Intensity  of  Toroidal  Coil    1 
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measured  values  are  "point"  values  only  insofar  as  the  field 
is  uniform  over  the  area  of  the  piok  up  ooil.   The  field  is 
not  uniform  since  it  falls  off  at  something  like  r**^  to  r~4 
depending  on  which  component  is  considered.   This  error  is 
certainly  not  the  limiting  error  of  the  measurements. 

The  sense  relation  of  the  various  components  could  not 
be  detected  with  the  pick  up  coil  and  voltmeter  alone.  However, 
using  the  oscilloscope  and  introducing  the  400  oycle  signal 
from  the  power  supply  as  a  sweep  synchronizing  signal  would 
resolve  the  180  ambiguity.   The  sense  can  also  be  easily 
determined  by  the  magnetometer  measurement  method  where  the 
vector  addition  to  the  earth's  field  is  immediately  apparent 
as  a  plus  or  minus  change • 

4.   Measurements  with  Field  Compensation 

As  a  possible  means  of  further  reducing  the  field  a 
"wand"  was  constructed  which  consisted  of  a  coaxial  cable, 
carrying  the  same  current  as  the  double  solenoid  coil,  on 
the  end  of  which  was  a  single  turn  loop.   The  loop  diameter 
was  eight  inches.  By  moving  this  loop  to  an  experimentally 
determined  position  between  the  pick  up  coil  and  the  double 
solenoid  it  was  possible  to  reduce  any  one  component  of  the 
field  to  an  "in  the  noise"  level.  However,  the  other  two 
components  would  crop  up  several  times  larger  than  before. 
Using  a  larger  wand  would  have  provided  a  more  uniform 
field  and  better  cancellation  could  have  been  achieved.   In 
fact  a  Helmholtz  pair  might  have  been  used  except  for 
limitations  of  space  in  an  airborne  system.   This  scheme, 
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of  course,  would  only  supply  "one  way"  canoellation.   The 

second  double  solenoid  placed  at  the  spot  of  the  pick  up 
coil  would  require  a  differently  configured  Helmholtz  pair 
to  compensate  for  its  external  field.   Such  a  configuration 
would  provide  more  mutual  coupling  so  that  transients  and 
power  supply  ripple  would  be  increased  problems.   However, 
better  field  reduotion  might  be  so  achieved  if  coaxial 
leads  were  used  throughout  and  if  great  care  were  taken  to 
avoid  any  uncompensated  short  leads.   The  above  method 
might  be  referred  to  as  determining  the  field  pattern  and 
then  "shimming"  the  field  by  means  of  coils  placed  external 
to  the  main  coil  system.   This  technique  could  be  applied 
to  any  ooil  configuration,  the  double  solenoid,  toroid  or 
flux  ball. 
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CHAPTER  VI 
SYSTEM  DESIGN 

1.  General  # 

The  dual  ooil  system  was  constructed  using  two 
Varian  V-4900  station  magnetometers  as  the  basic  "building 
blocks.   The  dual  system  differs  from  the  standard  system 
in  that  a  precession  signal  is  available  for  field  measure- 
ment almost  continuously.  For  a  complete  description  of 
system  design  reference  should  be  made  to  Yowell  (14). 
Figure  6-1  is  a  simplified  block  diagram  of  the  system, 

2.  Slow  Rate  Section 

Switch  S^  signifies  that  the  souroe  of  signal  is 
alternately  available  from  one  ooil  or  the  other.  After 
two  seconds  the  signal  has.  decayed  so  that  the  shift  to 
the  other  ooil  is  made.   The  signal  counter  gates  the 
500  KC  pulses  to  a  500  KC  counter  which  is  read  out  by 
an  analoging  circuit  every  4096  cycles  of  the  precession 
signal.   Switch  S2  is  closed  at  the  end  of  the  count  to 
store  the  analog  potential  on  capaoitor  C.   Sg  is  reopened 
after  about  60  milliseconds  so  that  the  next  count  may 
prooeed  without  ohanging  the  analog  potential  of  the  previous 
reading.  # 

3.  Fast  Rate  Section 

In  addition  to  the  4096  interval  count  a  512  interval 
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count   was   made   to   determine    the    field  strength  with  less 
aocuraoy  but  at  a  faster  rate.      The  advantage  of  having  a 
fast  rate   seotion  when  moving  rapidly  by  an  anomoly   is 
well   shown  by  figure  7  -   4D  and  E. 
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CHAPTER  VII 
SYSTEM  TEST 

1.  General 

The  system  was  placed  in  operation  with  the  coils 
separated  by  four  feet  as  shown  in  figure  7-1.   There  was 
a  fairly  high  gradient  of  the  earth!s  field  in  the  test 
area  due  to  the  presence  of  buildings.   Therefore,  it  was 

necessary  to  align  the  coils  along  a  constant  field  line 
so  that  both  coils  would  read  the  same  field  value. 
Actually  the  gradient  was  allowed  to  compensate  for  the 
residual  field  effect  between  the  two  coils.   Once  this  was 
done  the  tests  proceeded  without  changing  the  coil  positions. 

The  system  tests  were  made  using  kerosene  as  a  sample* 
due  to  the  insulation  problem  with  water.   Since  Tg  of 
kerosene  is  about  .7  seconds  the  slow  rate  count  interval 
was  changed  from  4096  to  a  2048  interval. 

2.  Signal  to  Noise  Ratio 

The  initial  signal  to  noise  ratio  of  the  system  was 
measured  as  70  to  one  (both  voltages  rms).   With  this 
condition  exieting  the  records  of  figure  7-2  were  taken. 
Analysis  of  these  show  that  the  fast  rate  has  an  rms 
error  of  .9  gamma  while  the  slow  rate  has  a  .31  rms  gamma 
error. 

The  sources  of  the  errors  are  the  least  count 
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error  (.1  gamma),  the  magnetic  noise  of  not  too  distant 
equipment,  thermal  noise*,  system  noise  and  analoging  errors. 
There  is  also  the  possibility  that  the  earth's  field  contains 
some  frequency  components  which  this  system  measures  and  of 
which  there  is  little  or  no  previous  knowledge.   As  an 
indication  of  this  possibility  during  the  period  of  this 
work  earth's  field  changes  were  recorded  at  the  rate  of 
10  gammas  per  minute  with  a  total  change  of  80  gammas  in 
eight  minutes.   This  rapid  change  correlated  with  sun 
spot  activity  and  observance  of  northern  lights. 

3.   Simulated  Anomoly  signal  tests 

In  order  to  simulate  an  aircraft  passage  over  an 
anomoly  two  tests  were  used.   One  was  driving  a  car  by  the 
system  at  two  distances,  30  feet  and  85  feet.   The  second  was 
to  rapidly  twirl  a  3  foot  iron  pipe  through  360  degrees  when 
at  a  distance  of  30  feet.   Lata  was  obtained  at  both  the 
slow  and  fast  rates.   The  records  show  a  short  dashed  line 
for  the  value  of  the  readings.   Every  other  line  is  from  the 
same  coil  at  the  2048  interval  rate.   Signals  from  the  same 
coil  are  indicated  by  the  marker  trace  at  the  bottom  of 
each  of  the  records.   For  the  fast  rate  every  other  four 
readings  is  from  the  same  coil.   Another  type  of  display 
is  given  in  the  recordings  and  this  is  an  apparently  contin- 
uous signal  obtained  by  feeding  the  discrete  signals  through 

a  low  pass  filter.   This  has  considerable  advantage  from  an 

operator's  point  of  view. 
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From  figure  7-3  it  is  seen  that  a  car  passing  at  30 
feet  creates  about  a  12  gamma  field  at  the  ooile.   Figure 
7-5  B  and  C  are  repeats  of  figure  7-3  C  and  D.   This  has 
been  done  to  demonstrate  the  effect  of  a  scale  shift.  This 
shows  that  the  system  does  not  saturate  but  will  shift 
scales  as  necessary  to  follow  the  oh&nging   field. 

One  very  interesting  aspect  of  figure  7-5  B  is  that 
it  clearly  shows  the  gradient  of  the  field  existing  between 
the  two  coils.   The  gradient  is  caused  by  the  fact  that  the 
oar  had  to  be  parked  at  a  distance  of  45  feet  from  the  coils 
between  rune  but  closer  to  one  coil  than  the  other.  All 
of  the  car  at  30  feet  signals  display  this  effect  except 
where  the  low  pass  filter  is  used.  The  gradient  here 
detected  is  about  .5  gammas  per  foot. 

Figure  7-4  B  and  C  shows  the  signal  from  a  car  passing 
at  85  feet  from  the  detection  coils.  This  is  about  a  one 
gamma  peak  signal  and  is  easily  disoernable. 

Figure  7-4  D  and  E  show  the  advantage  of  having  a  fast 

rate  count.  Figure  7-4  D  shows  a  filtered  2048  count  signal 
while  figure  7-4  E  shows  a  filtered  512  count  signal.   The 
height  of  the  peaks  is  at  least  doubled  by  the  fast  count 
system.  For  low  level  aerial  survey  work  the  fast  rate 
would  be  essential. 
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CHAPTER  VIII 
CONCLUSIONS 

The  fast  and  continuous  rate  magnetometer  has  been 

demonstrated  to  be  feasible.   The  one-tenth  gamma  accuracy 
initially  designed  for  was  not  achieved  primarily  because 

of  the  least  count  and  noise  errors.  The  least  count 

error  is  easily  oorreoted  by  multiplying  to  a  higher  frequency. 

Increasing  the  signal  to  noise  ratio  will  not  be  so  easy. 

To  increase  the  effective  signal  to  noise  ratio  the 
following  should  be  considered.  The  coil  design  can  be 
optimized  still  maintaining  the  requirements  for  low 
external  field.  This  inoludes  using  as  large  a  coil  as 
feasible  under  the  current  turn  off  restrictions.  The 
polarizing  current  may  be  increased  if  better  coil  design 
is  achieved  or  if  the  field  effect  between  the  two  coils 
is  compensated  for  by  other  means. 

This  compensation  might  be  done  by  determining  the 
effeot  of  position  modulation  (change  of  azimuth  of  the 
bird  carrying  the  coils)  on  the  readings  from  the  two 
ooils.   Qnoe  this  is  determined  it  could  be  fed  back  as 
a  compensating  .signal  to  make  the  two  coils  always  present 
the  same  reading  in  a  uniform  field.  The  compensation 
would  be  different  in  areas  where  the  earth's  field  vector 
has  changed  appreciably.   Such  a  compensation  system  would 
certainly  add  many  complexities  but  might  be  necessary  in 
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order  to  compensate  for  the  field  effect  of  the  carrying 
aircraft. 

A  simpler  means  of  achieving  this  compensation  would 
be  an  operator  balance  control.   This  would  provide  a  fine 

oontrol  of  the  analoging  potential  of  one  of  the  systems. 

The  operator  could  then  balance  the  readings  from  the  two 

coils  at  the  beginning  of  each  run. 

As  another  alternative  accept  the  display  showing  an 

effective  gradient  between  the  coils  due  to  the  external 
field  effect.   Where  absolute  field  measurements  are  desired 
the  data  reduction  process  could  make  the  necessary  corrections. 

There  is  a  possibility  of  finding  a  substance  which 
will  possess  a  short  T,  so  that  a  large  polarization  may  be 
obtained  in  the  same  time  of  polarizing  but  which  will  have 
a  longer  T2  so  the  signal  does  not  decay  so  rapidly.   Such 
a  short  T^  could  only  occur  if  the  lattice  energy  spectrum 
of  the  unknown  substance  contained  frequency  components 
corresponding  to  the  value  of  the  polarizing  field  as  related 
by   u)  -  ^>  H.   That  such  a  substance  may  exist  is  by  no 
means  a  certainty. 

Another  possibility  of  enhancing  the  signal  is  to 
synchronize  the  polarizing  current  turn  on  so  as  to  catch 
the  remaining  polarization  at  the  end  of  the  count  period. 
This  would  allow  the  polarization  to  build  up  to  a  higher 
value  for  a  given  polarize  time.   That  such  a  technique 
was  feasible  was  demonstrated  but  not  incorporated  in 
the  system. 
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The  coils  designed  and  used  demonstrated  primarily 
that  it  is  possible  to  operate  two  detection  systems  in 
close  proximity.   The  gradient  of  the  external  field  across 
the  signal  coil  does  not  appreciably  reduce  T2  and  readings 
are  therefore  possible. 

To  further  reduce  the  external  field  effect  the  coils 
oan  be  separated  a  greater  distance.  For  example,  separating 
the  double  solenoids  by  six  feet  would  reduce  the  external 
field  to  .1  gammas  per  ampere. 

To  improve  the  low  external  field  coil  design  the 
following  should  be  considered,  A  toroid  constructed  with 
a  large  overall  diameter  (12  inches  or  greater)  but  a  small 
torus  diameter  (one  inch)  would  greatly  roduoe  winding 
errors,  A  system  of  field  compensation  using  Helmholtz 
coils  could  be  used.  A  flux  ball  construction  could  be 
tried. 
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